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\l ! Abstract 

In recent years, the so-called polynomial moment problem, motivated 
by the classical Poincare center-focus problem, was thoroughly studied, 
' and the answers to the main questions have been found. The study of a 

, similar problem for rational functions is still at its very beginning. In this 

paper, we make certain progress in this direction; namely, we construct 
an example of a Laurent polynomial for which the solutions of the cor- 
, responding moment problem behave in a significantly more complicated 

way than it would be possible for a polynomial. 



1 Introduction 

The main result of this paper is a construction of a particular Laurent poly- 
nomial with certain unusual properties. This Laurent polynomial is a counter- 
■ example to an idea that, so far as the moment problem is concerned, rational 



functions would behave in the same way as polynomials. The main interest of 
the paper, besides the result itself, lies in a peculiar combination of methods 
I which involve certainly the complex functions theory but also group representa- 

■ tions, Galois theory, and the theory of Belyi functions and "dessins d'enfants" , 

while the motivation for the study comes from differential equations. 
- In addition to theoretical considerations our project involves computer cal- 

K> ' culations. It would be difficult to present here all the details. However, we tried 

^ . to supply an interested reader with sufficient number of indications in order for 

him or her to be able to reproduce our results. A less interested reader may 
omit certain parts of the text and just take our word for it. 
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About a decade ago, M. Briskin, J.-P. Frangoise, and Y. Yomdin in a series 
of papers [2]-[5] posed the following 

Polynomial moment problem. For a given complex polynomial P and dis- 
tinct complex numbers a, b, describe all polynomials Q such that 

b 

P'dQ ^ (1) 

for all integer i > 0. 

The polynomial moment problem is closely related to the center problem for 
the Abel differential equation in the complex domain, which in its turn may be 
considered as a simplified version of the classical Poincare center-focus problem 
for polynomial vector fields. The center problem for the Abel equation and the 
polynomial moment problem have been studied in many recent papers (see, e. g., 

m-i, m, m-m)- 

There is a natural sufficient condition for a polynomial Q to satisfy (H]). 
Namely, suppose that there exist polynomials P, Q, and W such that 

P = PoW, Q = QoW, and W{a) = W{b), (2) 

where the symbol o denotes a superposition of functions: /10/2 = /i(/2)- Then, 
after a change of variables z — > W{z) the integrals in ([1]) are transformed to the 
integrals 

/ P' dQ (3) 

Jw{a) 

and therefore vanish since the polynomials P* and Q are analytic functions in C 
and the integration path in ([3]) is closed. A solution of H]) for which ([2]) holds 
is called reducible. For "generic" collections P, a, b any solution of ^ turns 
out to be reducible. For instance, this is true if a and 6 are not critical points 
of P, see [I^, or if P is indecomposable, that is, if it cannot be represented 
as a superposition of two polynomials of degree greater than one, see |15| (in 
this case ^ reduces to the equalities P = W, Q = Q o P, and P{a) = P{b)). 
Nevertheless, as it was shown in jl4], if P{z) has several composition factors W 
such that W{a) — W{b) then the sum of the corresponding reducible solutions 
may be an irreducible one. 

It was conjectured in |16j that actually any solution of ([T]) can be represented 
as a sum of reducible ones. Recently this conjecture was proved in [T3]. The 
proof relies on two key components. The first one is a result of |T7] which states 
that Q satisfies ^ if and only if the superpositions of Q with branches P~^{z), 
1 < i < n, oi the algebraic function P^^(z) satisfy a certain system of linear 
equations 

n 

J2U^Q{P^^\^))=Q, fs,:eZ, 1 < 5 < fc, (4) 
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associated to the triple P, a, b in effective way. 

The second key componenet is related to the vector subspace Vp,a,6 C Q" 
spanned by the vectors 

(/s, £7(1)7 /s,(j(2), /s,(j(n)), 1 < S < fc, aeGp, 

where Gp is the monodromy group of P and /«, 1 < s < fe, are vectors from (|4]). 
By construction, the subspace Vp,a,b is invariant under the action of Gp, so the 
idea is to obtain a full description of such subspaces. In short, it was proved 
in [13] that if a transitive permutation group G < S„ contains a cycle of length n 
then the decomposition of Q" in irreducible components of the action of G 
depends only on the imprimitivity systems of G. Obviously, the monodromy 
group of a polynomial of degree n always contains a cycle of length n which 
corresponds to the loop around infinity. Furthermore, imprimitivity systems of 
Gp correspond to functional decompositions of P. Therefore, the structure of 
invariant subspaces of the permutation representation of Gp over Q depends 
only on the structure of functional decompositions of P, and a careful analysis 
of system dU and of the associated space Vp,a,6 eventually permits to prove 
that any solution of ([T]) is a sum of reducible solutions. Notice that using 
the decomposition theory of polynomials one can also show that actually any 
solution of ll]) may be represented as a sum of at most two reducible solutions, 
and describe these solutions in a very explicit form (see [10]). 

For example, in the simplest case of the problem corresponding to an inde- 
composable polynomial P the above strategy works as follows. The only invari- 
ant subspaces of the permutation representation of Gp on Q" in this case are 
the subspace Vi spanned by the vector (1,1,...,!), and its complement V-^^. 
Since system Q contains an equation whose coefficients are not all equal this 
implies that Vp^a,b = and therefore ([4]) yields that 

QiP,-\z)) ^ QiP,-\z)) = • • • = QiP-\z)) 

identically over z. On the other hand, such an equality is possible only if 
Q = Q o P for some Q G C[z]. Finally, P(a) — P{b) since otherwise after the 
change of variables z P{z) we would obtain that Q is orthogonal to all powers 
of z on [P(a), P(6)] in contradiction to the Weierstrass theorem. 

In the paper |19j the following generalization of the polynomial moment 
problem was investigated: for a given rational function F and a curve 7 C CP^, 
describe rational functions H such that 

F'dH ^ (5) 

for all i > 0. In particular, in [T9' another version of system ([4]) was constructed: 
its solutions, instead of the equality (O, guarantee only the rationality of the 
generating function f(t) — X^i^o '^i^* for the moments 

m, = / F'dH . (6) 
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On the other hand, it was shown that if the additional conditions H~^{oo} C 
F~-'^{(X)} and F{oo) — oo are satisfied, then the rationaUty of f{t) actually 
implies that f{t) = 0. 

The following modification of ([2]) is a natural sufficient condition imply- 
ing ([5]): there exist rational functions F, H, and W such that 

F ^ FoW, H = HoW , (7) 

the curve M^(7) is closed, and all the poles of the functions F, H lie "out- 
side" VF(7) (the term "outside" is written in quotation marks since it is defined 
also for self-intersecting curves). We will call such a solution of ([5]) geometri- 
cally reducible. Notice that if 7 is closed then geometrically reducible solutions 
always exist. Indeed, one may take 

W = F, H = HoF 

where H is any rational function with all its poles outside the curve F{-^). 
It is also shown in [THj that, similarly to the case of a polynomial P, for a 
generic rational function F (for example, for an F whose monodromy group 
is the full symmetric group) all solutions of ([5]) turn out to be geometrically 
reducible. However, for a non-generic F the situation becomes much more 
complicated in comparison with the polynomial moment problem, and some 
reasonable description of solutions of ([5]) seems (at least for the moment) to be 
unachievable. 

In this paper we will consider a particular case of problem ([5]) which is 
especially interesting in view of its connection with the classical version of the 
Poincare center-focus problem. Namely, we will consider the following 

Laurent polynomial moment problem. For a given Laurent polynomial L 
which is not a polynomial in z or in l/z, describe all Laurent polynomials Q 
such that 

[ L'dQ ^ (8) 

for all integer i > 0. 

In contrast to the polynomial moment problem, not any solution of the Lau- 
rent polynomial moment problem is a sum of geometrically reducible solutions. 
For example, as it was observed in 19J, if L{z) = L{z'^) for some d > 1, then 
the residue calculation shows that condition ([8]) is satisfied for any Laurent 
polynomial Q containing no terms of degrees which arc multiples of d. We will 
call such a solution of the Laurent polynomial moment problem algebraically 
reducible. Notice that, in distinction to geometrically reducible solutions which 
always exist, algebraically reducible solutions exist only if L is decomposable 
and has z'' as its right composition factor. One might think that any solution 
of the Laurent polynomial moment problem is a sum of geometrically and/or 
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algebraically reducible solutions but, as we will see below, this is not the case 
either, although it seems that for a "majority" of Laurent polynomials L this is 
the case. 

It is natural to start the investigation of the Laurent polynomial moment 
problem by the study of the particular case where L is indecomposable. At 
least, in this case there exist no algebraically reducible solutions. On the other 
hand, any geometrically reducible solution of ^ must have the form Q = QoL, 
where Q is a rational function whose poles lie outside the curve L(S^). However, 
since Q is a Laurent polynomial, it is easy to see that in this case Q is necessarily 
a polynomial. Furthermore, a sum of geometrically reducible solutions has the 
form 

and hence is itself a geometrically reducible solution. Thus, "expectable" (and 
therefore not very interesting) solutions of the Laurent polynomial moment 
problem for indecomposable L are of the form Q = Q{L), where Q is a poly- 
nomial. Any other solutions, when they exist, are of great interest since they 
show that the situation is more complicated than one might hope. 

Let L be an indecomposable Laurent polynomial of degree n, and let Gl be 
its monodromy group. We will always assume that L is proper, that is, it has 
poles both at zero and at infinity. In this case the group Gl contains a permuta- 
tion with two cycles: this permutation corresponds to the loop around infinity. 
Furthermore, it follows from Theorem 4.5 of [H] that if the only invariant sub- 
spaces of the permutation representation of Gl on Q" are Vi and Vf'-, then 
any solution of the Laurent polynomial moment problem for L is geometrically 
reducible. Therefore, if we want to find an example of a Laurent polynomial L 
for which there exist solutions which are not geometrically reducible, we may 
use the following strategy: 

• First, find a permutation group G of degree n such that G would contain 
a permutation with two cycles, and the permutation representation of G 
on Q" would have more than two invariant subspaces. 

• Then, realize G as the monodromy group of a Laurent polynomial L. 

• And, finally, prove somehow the existence of non-reducible solutions. 

This program was started in |19| . Namely, basing on Riemann's existence 
theorem it was shown that there exists a Laurent polynomial L of degree 10 such 
that its monodromy group is permutation isomorphic to the action of S5 on two- 
element subsets of the set of 5 points. The corresponding permutation action 
of S5 on Q^*^ has more than two invariant subspaces. Furthermore, proceeding 
from a general algebraic result of Girstmair [TTj about linear relations between 
roots of algebraic equations it was shown that there exists a rational function Q 
which is not a rational function in L such that the generating function for the 



5 



sequence of the moments 



SI 

is rational (see Sec. 8.3 of [H]). However, the methods of [Hj do not permit to 
find L 01 Q exphcitly and tell us nothing about the structure of solutions of ([5]) . 

In this paper we provide a detailed analysis of the above example with the 
emphasis on the two following questions of a general nature: 

• First, how to construct a Laurent polynomial L starting from its mon- 
odromy group Gl? 

• Second, how to describe solutions of ^ which are not geometrically re- 
ducible ? 

We answer all these questions for the particular Laurent polynomial L given 
below. Actually, we believe that our methods can be used in a more general 
situation too and can serve as a "case study" for further research concerning 
the Laurent polynomial moment problem. 

The main result of this paper is an actual calculation of an indecomposable 
Laurent polynomial L such that the corresponding moment problem has non- 
reducible solutions, and a complete description of these solutions. Namely, we 
show that for 

_ K{z-lf{z-af{z-b) 
L , (9) 

where 

115/- 3173 

K = \ V5, a = h - V 5, b = VS , 

216 216 ' 2 2' 22' 

there exist Laurent polynomials Qq = 1, Qi,Q2,Q3,Qi (we compute them 
explicitly in Sec. |4|) such that the following statement holds: 

Theorem 1.1. A Laurent polynomial Q is orthogonal to all powers of L on 
if and only if Q can be represented in the form 

4 

Q = Y,{R,oL)-Q, 

for some polynomials Rq, Ri, R2, R3, R4. 

In other words, solutions of the moment problem for L form a 5-dimensional 
module over the ring of polynomials in L (while in a generic case such a module 
is one-dimensional and is therefore composed of polynomials in L and of nothing 
else). The choice of the basis Qj is not unique, but once a basis is chosen the 
above representation of Q becomes unique. 
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The paper is organized as follows. In Sec. 2 we give a detailed description of 
the permutation action of S5 on Q"'^". In Sec. 3 we compute explicitly a Laurent 
polynomial L whose monodromy group is permutation equivalent to this action. 
Finally, in Sec. 4 we determine the above mentioned Laurent polynomials Qj 
and prove Theorem ll.il 
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2 Permutation representation of S5 on Q^^ with 
more than two invariant subspaces 

Consider the complete graph = [V^ E) having the vertex set V = {1, 2, 3, 4, 5} 
and the edge set E consisting of all the subsets of V of size 2. The symmetric 
group S5 acts on V and therefore also on i?, and we thus obtain a transitive 
action of S5 of degree 10. Moreover, the homomorphism S5 Sio is obviously 
injective. Let us identify the canonical basis 

ei = (1,0,0,0,0,0,0,0,0,0), 
62 = (0,1,0,0,0,0,0,0,0,0), 

ew = (0,0,0,0,0,0,0,0,0,1) 

of the space Q^*^ with the set E. This identification may in principle be arbitrary 
but we have chosen the one which is more "readable" , see Fig. [TJ the first five 
vectors are associated, in a cyclic way, to the sides of the pentagon, while the 
last five vectors are associated in the similar way to the sides of the inside 
pentagram. 

Associating to each element of S5 a lOx 10 permutation matrix corresponding 
to the action of this element on E we obtain a permutation representation of S5 
on Q^°. Any permutation representation of any finite group always has at least 
two invariant subspaces: the subspace Ui of dimension 1 spanned by the vector 
1 = (1,1,...,1), and its orthogonal complement t7„-i — of dimension n—1 
containing the vectors (xi, X2, . . . , Xn) having J27=i — While the space Ui 
is obviously irreducible, the space Un-i may be, or may not be irreducible. We 
will show that in our case it is reducible. 

One of the ways to construct invariant subspaces in our example is to consider 
subsets of edges which are sent to one another by the action of S5 on the vertices. 
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Figure 1: The correspondence between the edges of the complete graph and the 
basis in O^". 



Let us take the fans Fi C E, i = 1, . . . ,5, where Fi is the set of edges of 
incident to the vertex i, see Fig. [51 




Figure 2: The fans Fi, that is, the sets of edges incident to the vertex i = 1, . . . ,5. 
Obviously, any permutation of the vertices sends fans to fans. Therefore, 
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the vectors Vi = J2ueF- O'"' more concretely, 

Vl = (1,0,0,0,1,1,0,0,1,0), 

V2 - (1,1,0,0,0,0,1,0,0,1), 

n = (0,1,1,0,0,1,0,1,0,0), 

V4 = (0,0,1,1,0,0,1,0,1,0), 

V5 = (0,0,0,1,1,0,0,1,0,1) 

(the first five and the last five components of these vectors move cyclically) span 
an invariant subspace F C Q^°. It is easy to verify that F is 5-dimensional. 
Since every edge is contained in exactly two fans we have = {2,2, . . . , 2) 

and therefore F contains Ui as its subspace. The orthogonal complement of Ui 
in F is a 4-dimensional invariant subspace U4 C Q^". The vectors 

V2~vi ^ (0,1,0,0,-1,-1,1,0,-1,1), 

V3~vi = (-1,1,1,0,-1,0,0,1,-1,0), 

Vi-vi = (-1,0,1,1,-1,-1,1,0,0,0), 

V5-V1 = (-1,0,0,1,0,-1,0,1,-1,1), 

each having equal number of ones and minus ones, are orthogonal to the vector 1. 
They are linearly independent, and therefore they span U4. 

Another collection of subsets of E which is stable under the action of S5 
is the set of Hamiltonian cycles H C E, that is, cycles that visit each vertex 
exactly once. A Hamiltonian cycle in can be described by a 5-cycle c € S5 
which indicates in which order the vertices are visited; note that describes 
the same Hamiltonian cycle since our graph is undirected. The complement 
H = E \ H is also a Hamiltonian cycle which corresponds to the permutation 

(or to its inverse c~^). There are 24 cyclic permutations in S5; they give rise 
to 12 Hamiltonian cycles in Kc, which form 6 pairs of mutually complementary 
cycles: see Fig. [H 

The vectors Wk = Y^ueH^ " Y.ue'Hk ^" ^o^'s concretely, 

wi = (1,-1,1,-1,1,-1,1,1,-1,-1) 

W2 = (1,1,-1,1,-1,-1,-1,1,1,-1) 

W3 = (-1,1,1,-1,1,-1,-1,-1,1,1) 

Wi = (1,-1,1,1,-1,1,-1,-1,-1,1) 

W5 = (-1,1,-1,1,1,1,1,-1,-1,-1) 

we = (-1,-1,-1,-1,-1,1,1,1,1,1) 

(once again the first five and the last five components move cyclically) span 
an invariant subspace. Every edge of belongs to 3 "positive" Hamiltonian 
cycles and to 3 "negative" ones; therefore, X]i=i '^j = 0. It is easy to verify that 
the space C/5 spanned by these 6 vectors is in fact 5-dimensional. For every fan 
Fi and for every pair {Hj, Hj), exactly two edges of Fi belong to Hj, while the 



9 



1 1 1 




i?4 He 

Figure 3: Pairs of Hamiltonian cycles in Kr,: the "positive" cycles Hi arc drawn in 
bold lines, while their complements, the "negative" cycles Hi, are shown in thin lines. 

other two belong to Hj. Therefore, Vi _L Wj for all so U5 _L F where, as 
before, F ^Ui® U4,. 

Thus, we get a decomposition of into three invariant subspaces: 

= Ui (B U4 (B U5. We did not prove that the subspaces U4 and U5 are 
irreducible. The proof goes by some routine verification using the character ta- 
ble of S5. We omit the details since for our goal this fact is irrelevant: the only 
thing we wanted to show was the reducibility of the orthogonal complement 
- = Ug, and this statement is proved since we have shown that Ug = 1/4® U5. 

We finish this section by specifying how certain elements of S5 act on the 
labels of the 10 edges. By construction, the permutation / = (1,2,3,4,5) G S5 
acts as 

<^ = (1,2,3,4,5)(6,7,8,9, 10). 
Taking a simple transposition, for example, a = (2, 5) G S5, we get 

a=(l,5)(2,8)(4, 7). 

Indeed, all the edges having both ends different from 2 and 5, remain fixed, as 
well as the edge {2, 5} itself, while the 6 edges having exactly one end equal to 
2 or to 5 split into 3 pairs. Finally, taking s = (1,2)(3,5,4) e S5 we obtain 

<7=(2,5,7,6, 10,9)(3,8,4). 
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Note that — (1, 2); conjugating this element by / we get ah the transpositions 
(i, ? + 1) of adjacent elements. Therefore, the elements and /, and hence also 
s and /, generate the whole group S5. Since 

aaip = 1 

and the homomorphism S5 Sio is injective, this implies that the group 
{a,a,^p) is generated by a and a and is isomorphic to S5. The action of 
(tj, a, ip) = S5 on the 10 edges is primitive; indeed, we could only have 2 blocks 
of 5 elements each, or 5 blocks of 2 elements each, but the presence of a cy- 
cle of order 6 is incompatible with the first possibility while the presence of a 
single fixed point is incompatible with the second one. The action is obviously 
transitive. 

3 Realization of the degree- 10 action of S5 as the 
monodromy group of a Laurent polynomial 

During all this section, we systematically use various methods and results of 
the theory of "dessins d'enfants". We will try to be concise but clear. For all 
missing details the reader may address the book 12J (Chapters 1 and 2). 

3.1 Belyi functions and "dessins d'enfants" 

Rational functions from CP^ to CP^ (and, more generally, meromorphic func- 
tions from a Riemann surface X to CP^), unramified outside 0, 1, and oo, are 
called Belyi functions. They have many remarkable properties. In particular, 
any such function F{x) may be "encoded" in the form of a bicolored map Mp 
drawn on the sphere (resp., on the surface X). Namely, let us color the points 
and 1 in black and white respectively, draw the segment [0, 1], and define Mp 
as the preimage Mp = i^~^([0, 1]) of the segment [0, 1] with respect to the func- 
tion F{x) : CP""^ CP^. By definition, black (resp., white) vertices of Mp are 
preimages of the point (resp., of the point 1) and edges of Mp are preimages 
of the segment [0, 1]. 

The segment [0, 1] may itself be considered as a bicolored map having two 
vertices of degree 1 and a face of degree 2 containing infinity. Clearly, Mp 
has n = deg F edges, and the degree of a vertex x of Mp coincides with the 
multiplicity of x with respect to F. Furthermore, each face of Mp contains a 
pole of F, and twice the multiplicity of this pole coincides with the degree of 
the corresponding face. The map Mp permits to reconstruct the monodromy 
group Gp of F. Indeed, let go, 9i be generators of Gp corresponding to the 
loops around and 1. Taking a base point of the covering somewhere inside the 
segment [0, 1] we may consider that the permutations and gi act not on the 
preimages of the base point but on the preimages of [0, 1], that is, on the edges 
of Mp. The permutation (resp., 171) sends an edge e to the next one in the 
counterclockwise direction around the black (resp., white) vertex adjacent to e. 
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Notice that if goo is the element Gp corresponding to the loop around oo, 
then gogigoc = 1- 

For example, assuming that a Belyi function F corresponds to the map shown 
in Fig. SI we may conclude that F is of degree 10 (since there are 10 edges), has 
two poles, both of order 5 (since there are two faces, both of degree 10), and 
that the corresponding permutations go,gi,goo coincide with the permutations 
cr, a, and ip defined at the end of the previous section. 




Figure 4: Realization of Ss acting on 10 edges of a bicolored plane map. 

Riemann's existence theorem implies that for any bicolored plane map there 
exists a Belyi function F{x) which is unique up to a composition with x l—^ fj.{x) 
where /i(x) is a linear fractional transformation. In particular, since for the map 
shown in Fig.|3]the permutations go, gi, g^c coincide with ct, a, ip, this pictures 
"proves" that there exists a rational function F{x) whose monodromy group is 
permutation equivalent to the action of S5 on 10 points discussed above. Our 
next goal is to find this function explicitly. 

3.2 A system of equations for the coefficients of Belyi 
function, and its solutions 

In the rest of this section we will compute a Belyi function which produces a 
map isomorphic to that of Fig. [4] as a preimage of the segment [0, 1]. A reader 
not interested in the details of the computation may just take our word for it 
that the resulting function is the one given in ([§]), and pass directly to Sec. 4. 
We will provide not all the details but only a minimum allowing the reader to 
reproduce our results. 

The black vertices of the map are the preimages of 0, or, in other words, 
they are roots of the rational function F we are looking for. Furthermore, the 
vertex of degree 6 is a root of multiplicity 6, the vertex of degree 3 is a triple 
root, and the vertex of degree 1 is a simple root. The freedom of choosing a 
linear fractional transformation fj.{x) allows us to put these three points to any 
three chosen positions. Let us put, for example, the vertex of degree 6 to a; = 0, 
the vertex of degree 3, to a; = 1, and the vertex of degree 1, to x = — 1. Then, 
the numerator of F will take the form x^{x — l)^{x + 1). 
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The permutation corresponds to the monodromy above oo, and it has two 
cycles of length 5. Therefore, the function in question must have two poles of 
degree 5, one pole inside each face of the map. Suppose these poles to be the 
roots of a quadratic polynomial -[-ax + h. Then, the Belyi function in question 
takes the form 

(x'^ + ax + o)° 

where a, b are constants that remain to be determined. 

Here the reader may be surprised. We are looking not for an arbitrary Belyi 
function but for a Laurent polynomial, aren't we? Then, would it not be a 
better idea to use the same liberty of choice of three parameters and to put one 
of the poles to a; = 0, and the other one, to x = oo? The answer is no: such 
a choice would not be a good idea - at least at this stage of the computation. 
The reason is related to Galois theory and will be explained later, in Sec. 13.41 

The white vertices of our map are the preimages of 1, or, in other words, the 
roots of the function F{x) — 1. There are three white vertices of degree 2; they 
correspond to double roots of F{x) — 1. Computing the derivative of F we get 



x^{x — p{x) 
(x2 + ax + h)^ 



F'{x) = K 
where 

p{x) = (5a + 2)x^ + (2a + 106 + 4) x^ - (a - 26) x - 66 . 

It becomes clear that p{x) is the cubic polynomial whose roots are the three 
white vertices of degree 2, so the numerator of F{x) — 1 must have p(x)^ as a 
factor. Note also that the leading coefficient of this numerator is K — 1. Thus, 
we can now write down the hypothetical form of F{x) — 1 which we temporarily 
denote by Fl(x): 

^ K-1 p{xfq{x) 
(5a + 2)2 (a;2 + ax + 6)5 

where q{x) is yet unknown polynomial of degree 4, with the leading coefficient 1, 
whose roots are the four white vertices of degree 1. Denote 

q{x) — x"^ + cx^ + dx^ + ex + f 

and compute the derivative of H{x). 

The results of the subsequent computations become more and more cum- 
bersome. Their main steps go as follows. First of all, H{x) is nothing else but 
another representation of F{x) — 1, so we must get in the end F'{x) = H'{x). 
Therefore, after having computed H'{x) we ask Maple to factor the difference 
F'{x) — H'{x), and we get an expression 

F'ix)~H'{x) = Const ?'(^) ^(^) 



2 + ax + 6)6 

where r(x) is a (very huge) polynomial of degree 7. The final action to do is to 
equate r(x) to zero: this means that we extract its coefficients and equate all 
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of them to zero. This gives us a system of algebraic equations on the unknown 
parameters K, a, b, c, d, e, f. 

The solution of the system thus obtained using the Maple-7 package takes 
14 seconds. It takes significantly more time to enter all the involved formulas 
and operations. And it takes even more time to find our way among the solutions 
since they are many and varied. 

3.3 Finding our way among the solutions 

3.3.1 Maps v^rith the same set of vertex and face degrees 

If we analyse carefully the above procedure of constructing a system of equa- 
tions, we will see that the only information we have used about the map of Fig. [5] 
is the set of degrees of the black vertices, the white vertices, and the faces of 
this map. However, there exist not one but 7 maps having the degree partition 
of the black vertices equal to (6, 3, 1) = 6^3^ 1^ h 10, that of white vertices equal 
to (2, 2, 2, 1, 1, 1, 1) = 2^1-^ h 10, and that of the faces equal to (5, 5) = 5^ h 10. 
These maps are shown in Fig. [5l Therefore, the above computation must pro- 
duce Belyi functions for all of them. 




Figure 5: All the seven bicolored maps with the degree partition of the black vertices 
being 6^3^!^, that of the white vertices being 2'^1'*, and that of the faces being 5^. 
Monodromy groups are also indicated. 



The picture convinces us that these 7 maps do exist. In order to prove 
that there are no others we may compute the number of triples of permutations 
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{9o,9i,9oo) of degree 10 having the same cycle structure as {a,a,(p) and sat- 
isfying the equahty (7o5i5oo = 1- For this end, we may use, for example, the 
following formula due to Frobenius: 

Proposition 3.1. Let Ci, C2, . . . ,C'k be conjugacy classes in a finite group G. 
Then the number J\f{G;Ci,C2, ■■■ ,C'k) of k-tuples (xi, X2, ■ ■ ■ , Xk) of elements 
Xi Cz G such that each Xi G Gi and X1X2 ■ ■ ■ Xk = 1, is equal to 

N{G,GuG2,...,Cu) - ^ X. (dimx)^-2 ' 

where the sum is taken over the set of all irreducible characters of the group G. 

Applying this formula to the group G — Sio, k — and the conjugacy 
classes Ci, C2, C3 determined by the cycle structures 6^3^1^, 2^1^, and 5^, 
respectively, and computing the irreducible characters of Sio using the Maple 
package combinat, we get 

M{G;Gi,G2,Cs.) = 25 401600 = 7-10!. 

None of the maps shown in Fig. [5] has a non-trivial orientation preserving 
automorphism; therefore, each of them admits 10! different labelings. 

It is useful to determine monodromy groups of the functions corresponding 
to the above maps. For the map in the upper left corner we know already that, 
by construction, it is isomorphic to S5. For the 5 maps shown in the lower part 
of the figure, the order of the group (which can be calculated by the Maple 
package group, function grouporder) is equal to 10!, and therefore the group is 
Sio itself. Finally, for the map in the upper right corner, using the same Maple 
package, or GAP, or the catalogue [9], we may establish that it is isomorphic 
to Se. 



3.3.2 Galois action on maps and finding F{x) 

We find the coefficients of the Belyi functions by solving a system of algebraic 
equations. Therefore, there is no wonder that these coefficients are algebraic 
numbers. The group of automorphisms of the field Q of algebraic numbers is 
called the absolute Galois group and is denoted by F = Gal(Q|Q). An element 
of the group F, acting simultaneously on all the coefficients of a given Belyi 
function, transforms it into another Belyi function which may correspond to 
another map. 

Thus, bicolored maps split into the orbits of the above Galois action. The 
set of degrees of black and white vertices and faces is an invariant of this action; 
therefore, all the orbits are finite. Another invariant is the monodromy group. 
Looking once again at Fig. [5] we see that the set of 7 maps represented there 
splits into at least three Galois orbits: two orbits contain each a single element, 
while the set of the remaining 5 elements may constitute one orbit or further 
split into two or more orbits. The general theory suggests that for the singletons 
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the coefficients of the corresponding Belyi functions must be rational numbers. 
And indeed, among our solutions we find two such functions: 



Fi{x) 



50000 + 1) 



27 



and 



F2{x) = 337500- 



+ 4a; - 



x'^ix- + 1) 

(lla;2 + 4a;_i6)5 



At this stage we simply ask Maple to draw the F-preimages of the segment [0, 1] 
and find out that the function we are looking for is Fi: just compare Fig. [S] with 
Fig. m It is pictures like that in Fig. [6l obtained as Belyi preimages of the 
segment [0, 1], which are usually called dessins d'enfants. 



0.4 
0.2 


-0.2 
-0.4 



-1 



-0.5 







0.5 



1 



1.5 



Figure 6: A Maple plot of the "dessin d'enfant" corresponding to the Belyi func- 
tion Fi. Black vertices are marked by little squares. 



The five remaining maps constitute an orbit of degree 5 defined over the 
splitting field of the polynomial 

Q{t) = 85237 - 95206 + 48850 - 7456 f + 1606 i - 226 . 

This means that the coefficients of a Belyi function F(x) are expressed in terms 
of (more exactly, as polynomials of degree < 4 in) a root of this polynomial. 
Taking one by one five roots we obtain five different Belyi functions which 
correspond to the five maps with the monodromy group Sio shown in the lower 
part of Fig. [HI 

Notice that besides the solutions mentioned above, our system of algebraic 
equations produces a bunch of the so-called "parasitic solutions" representing 
various kinds of degeneracies. Some of them are easy to eliminate, others are 
not. For example, in one of the solutions we get a = 0, 6 = 0, which means that 
the denominator of F is x^^ , while its numerator contains x^ . This solution does 
correspond to a Belyi function, but of degree 4 instead of 10. Another easy case 
is K = 1, a = —2/5, which leads to a division of zero by zero in the constant 
factor (K — l)/(5a + 2)^ of the function H in Sec. 13.21 More difficult cases of 
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degeneracies also exist but we will not go here into further details, as well as into 
many other subtleties proper to any experimental work. The questions already 
discussed show quite well why the computation of Belyi functions remains a 
handicraft instead of being an industry. 

3.4 Prom a rational function to a Laurent polynomial 



Now we may return to the question asked in Sec. l3.2l and explain why we decided 
to compute a "generic" Belyi function instead of looking from the very beginning 
for a Laurent polynomial. 

We see that, while Fi is defined over Q, its two poles are not: they are 
roots of the quadratic polynomial x'^ + Ax — 1; concretely, they are equal to 
— 2 ± VS- Any linear fractional transformation of the variable x sending one 
of theses poles to and the other one to oo would inevitably add \/5 to the 
field to which belong the coefficients of Belyi functions. Thus, the functions 
defined over Q would become defined over Q(\/5), the orbit of degree 5 would 
become one of degree 10 (with each of the five maps being represented twice), 
parasitic solutions would also become more cumbersome (and their parasitic 
nature would be more difficult to detect), and so on. And without doubt Maple 
would have a much harder work to solve the corresponding more complicated 
system of algebraic equations. 

But from now on, after making all the above computations with the simplest 
possible fields, we can easily transform Fi into a Laurent polynomial. The 
transformation 

_ (2 - V5)x - 1 
^ ^ (2 + V5)x - 1 

sends the pole —2 — \/5 to and — 2 + Vs to oo, and also sends to 1. 
Substituting into -^1(2;) its inverse 

_ z-l 

(2 + V5)z - (2 - V5) 

we obtain the Laurent polynomial of Theorem 11.11 

iz-inz-afiz-b) 
L{z) = K 

where 

11 + 5V5 -3 + \/5 , 7-3^/5 

K = , a — , — . 

216 2 ' 2 

4 Proof of the main theorem 

We are looking for Laurent polynomials Qj, < j < 4, of the form 

i 

k=-j 
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(we set Qo = 1) satisfying the equation ([5]). However, it is clear that we may 
multiply Qj by a constant, and also add to Qj an arbitrary linear combination of 
Qi for i < j, and this gives us another solution having the same form. Therefore, 
in order to achieve uniqueness, we impose on Qj the following three conditions: 

1. The coefficient s_j is equal to one. 

2. For i = — j + 1, . . . , the coefficients Si are equal to zero. 

L'dQj =0 for aU l<i<j. 

The Laurent polynomial Qj has 2j + 1 coefficients; the first two conditions 
fix J + 1 of them, while the third condition provides us with j additional lin- 
ear equations on coefficients. In order to ensure that the integrals in the third 
condition vanish, according to the Cauchy theorem, we must calculate the co- 
efficients preceding z^^ in • Q'j, 1 < j < 4, 1 < i < j, and set them to zero. 
The existence and uniqueness of solutions will be explained later (in Step 3 of 
the proof). The results of the calculation are collected below. 

Qo - 1, 

_ z' + l 
Wl — , 



Q2 — 5 



Q3 



^47 I 2 
\ 2 ^ 2 



-(9 + 4 v^) -f (130 + 58 VE) z' - (630 + 282 ^5) -H (910 406 ^5) + 1 



Now, everything is ready in order to prove the main theorem. The proof is 
divided into several steps. 

Step 1. First of all, we must check that the Laurent polynomials Qj , 1 < j < 4, 
satisfy the equalities 

L'dQ = (11) 

for all i > (for Qo it is obvious). For this purpose we may use Theorem 7.1 
of [19] and verify this equality only for a finite number of i, namely, for 

1 < « < (A^- 1) -dcgQ + l, (12) 
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where N is the size of the orbit of the vector 

(1,1,1,1,1,-1,-1,-1,-1,-1) 

under the action of the monodromy group of L. In our case, degQj — 2j, 
1 < J < 4, and N = 12; therefore, the maximal value of the right hand side 
of (fT^ is equal to 89. The verification for all the four polynomials Qj takes less 
than one minute of work of Maple-11. 

Step 2. Observe that if a Laurent polynomial Q is a solution of pT|) then 
for any polynomial R the Laurent polynomial Q — R{L) ■ Q is also a solution 
of (fTTjl . Indeed, it is enough to prove it for i? = z'^, fc > 1. We have: 

f Vd{L^Q) = I V+^dQ+ ( VQdL^. 

<J Si <J Si ^ Si 

The first integral in the right-hand side of this equality vanishes by (HJ). On 
the other hand, for the second integral we have: 

/ L'QdL^ = k[ L'+'^-iQdL = / QdL^+' = ^/ L^+'dQ, 

Js^ Js^ t + k Jsi « + fc Jsi 

and therefore this integral also vanishes. 

Step 3. The final ingredient we need is Theorem 6.7 of 19 which states that 
if the leading degree of a Laurent polynomial L is a prime number p (in our 
case p — 5), and if Q is a polynomial (that is, a common one, not a Laurent 
polynomial) such that pT|) holds, then either L{z) = Li{zP) for some Laurent 
polynomial Li while Q is a linear combination of the monomials with i not 
being multiples of p, or Q is a constant. Since the Laurent polynomial L we 
are working with is not of the form L{z) = Li{z^) this result implies that a 
polynomial Q cannot satisfy (|lip unless Q is a constant. 

This fact also explains the uniqueness of Qj. Indeed, if o'iP and Q^-^-* are 
two solutions of the equations imposed on Qj at the beginning of this section, 

then their difference Q^^'' — Q^'p is also a solution of the Laurent polynomial 
moment problem. But this difference is a polynomial (since the terms z^^ in 
Qj^'' and Q''^ cancel) and therefore must reduce to its constant term; but the 
constant term of this polynomial is equal to zero. 

The uniqueness of the solution implies the non-degeneracy of the matrix of 
the system, and the non-degeneracy, in its turn, implies existence. 

Step 4. Now, let us suppose that Q is a Laurent polynomial satisfying (jlip 
and m{Q) < is the minimal degree of a monomial in Q. Let m{Q) — — 5fco— joi 
where < jo < 4 and fco > 0. Then for any cq G C the Laurent polynomial 
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is a solution ([TT|) . Furthermore, choosing an appropriate cq we can assume that 
m (Q(i)) > m{Q) (here we use the fact that the coefficient s-j in Qj is not zero). 
Now, if m(Q(^)) = — 5fci — ji, where < ji < 4 and ki > 0, then, setting 

g(2) ^ g(l) _ ^fei . g^.^ 

for an appropriate ci we obtain a solution (fTTjl with m{Q'^^^) > m((5'^^)). Con- 
tinuing in this way we will eventually arrive to a solution Q^*"' of (|lip for which 
m (QW) > 0. In view of the result cited in Step 3 such a solution should be a 
constant c e C. Therefore, 

r-l 4 
i=0 j=0 

for some polynomials Rq, Ri, R2, R3, Ri- The theorem is proved. 

Final remarks. In general, it is not known if the reducibility of the action of 
the monodromy group Gl of a Laurent polynomial L of degree n on the space Q" 
always implies a non-trivial structure of solutions of the corresponding moment 
problem. The only facts which follow from the general theory are as follows: 

• The reducibility of the above action implies the existence of a rational 
function Q, which is not a rational function in L, such that the generating 
function for the sequence of the moments 

^ VdQ, i> 0, 

is rational (see Sec. 8.3 of [ilO ). 

• If the above function Q turns out to be a Laurent polynomial, then the ra- 
tionality of the generating function implies its vanishing (see Theorem 3.4 

of US]). 

It would be interesting to understand in a more profound way what is the 
underlying mechanism which relates the structure of solutions of the moment 
problem for L with the structure of the representation of G^. 
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